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Limits of the uni-directional pulse propagation approximation 
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I apply the method of characteristics to both bi-directional and uni-directional pulse propagation 
in dispersionless media containing nonlinearity of arbitrary order. The differing analytic predictions 
for the shocking distance quantify the effects of the uni-directional approximation used in many 
pulse propagation models. Results from numerical simulations support the theoretical predictions, 
and reveal the nature of the coupling between forward and backward waves. 



Published as J. Opt. Soc. Am. B 24, 2363 (2007). This 
version contains an additional appendix containing the 
MOC predictions for a time-propagated theory. 



I. INTRODUCTION 

Most approaches to optical pulse propagation rely on 
an approximation where the fields only propagate for- 
wards. Even the recently derived extensions of typical 
propagation methods used in nonlinear optics (e.g. 
assume a complete decoupling between oppositely propa- 
gating fields to optimise the calculation. Moreover, those 




111) or 
), are of- 



based directly on Maxwell's equations 
the second order wave equation (e.g. 
ten simplified to work in the forward-only limit, where 
backward propagating fields are set to zero. This is de- 
spite directional decompositions of Maxwell's equations 
(e.g. [3, 9]) indicating that nonlinearity inevitably cou- 
ples the forward and backward waves together - and even 
creates a backward field if one is not present. Usually 
we assume that a forward wave will not generate a sig- 
nificant backward wave via the nonlinearity because the 
backward component is very poorly phase matched 1 . In 
contrast, deliberately trying to phase match the back- 
ward wave was suggested in the 1960's [l(|, and some 
progress has been made in materials in recent years 
(e.g. [H E EE El ) ; however such attempts are not the 
focus of this paper. 

Here I use the phenomenon of carrier wave shock- 
ing [l5| as a tool to probe the fundamental limits of a 
uni-directional model under the influence of the intrin- 
sically bi-directional nonlinear coupling. Carrier wave 
shocks are discontinuities in the electric (and magnetic 
field) profiles, and develop as the nonlinear effects dis- 
tort the waveform as it propagates. Assuming an initially 
uni-directional field, I compare analytic predictions for 
the shocking distance from uni- and bi-directional the- 
ories. These are based on the Method of Characteris- 
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1 If the forward field has a wavevector fco and evolves as 
exp(+ikoz), the generated backward component will evolve 
as exp(— ik^z). This gives a very rapid relative oscillation 
exp(— 2ikoz), which will quickly average to zero. 



tics (MOC) [l6(, and show a clear difference between the 
predicted shocking distances in the two cases. The (ex- 
act) bi-directional theory is based on the second order 
wave equation (as in 15|, 17|), whereas the (approximate) 
uni-directional theory is derived from the G ± directional 
fields 3 description. I support the theoretical results 
with pseudospectral spatial domain (PSSD) simulations 
H for both second order (x^) and third order (x^) 
nonlinearities. 

The results in this paper act as a bound on the valid- 
ity of one-dimensional optical propagation models using 
a uni-directional approximation. Since linear dispersion 
and finite nonlinear response times will typically dimin- 
ish any generation of a backward wave, it is clear that 
any model which can be assumed uni-directional on the 
basis of this paper will be more so in practise. These 
results do not tell us whether the uni-directional approx- 
imation would be more or less robust for models incorpo- 
rating transverse effects such as self-focussing, but they 
at least establishes a point of reference, valid for beams 
with a weak spatial variation. Note also that I consider 
only propagation within bulk media, since surfaces or in- 
terfaces can cause reflections, and clearly require a bi- 
directional model. 

Section |TT] briefly describes the MOC, and presents a 
prediction for the shocking distance in the bi-directional 
case due to simple nonlinearities of arbitrary order; sec- 
tion IIIII follows with analogous calculations using an ex- 
plictly uni-directional wave equation. Next, section IIVI 
discusses the analytic and numerical results, section fVl 
considers the role of the backward field, and then section 
PVTl presents some conclusions. 



II. BI-DIRECTIONAL MODEL 

Analytic formulae for the shocking distance in mate- 
rials with instantaneous response have been calculated 
for both third-order nonlinearities (by Rosen [HI) and a 
simple second order case (by Radnor [l8]). Both calcula- 
tions used the MOC to predict the formation of a value 
discontinuity in the field at certain points within the elec- 
tric field profile of a pulse or wave. Here I generalise 
the treatment to allow for an instantaneous perturbative 
nonlinearity of arbitrary order, following the calculation 
of Kinsler et al. E3l- Note that the calculation below 
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can easily be generalised to include a sum of nonlinear 
terms, if so desired [13]. 
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FIG. 1: The progressive distortion of an initially sinusoidal 
wave profile as the shocking distance is approached for (a) x 
and (b) \ nonlinearities in a dispersionless medium. In the 
spectral domain, the distortion corresponds to a build-up of 
significant quantities of higher order harmonic content. 

In this paper, I consider nonlinearities of the simplest, 
rather than the most general form, because then solutions 
using the MOC can be found. Ignoring the tensor nature 
of the coefficients, the displacement field is 



D 



E + x {1) E + X (m) E r ' 



(1) 



The usual second order wave equation for E in this case 
is 



, d 2 E 
dz 2 



, d 2 E 

1 dt 2 



■ X 



('") 



d 2 E" 



(2) 



where e r = 1 



X 



(1) 



no is the (relative) dielectric 
constant and no the linear refractive index. 

The MOC treats each point on the waveform E(t) sep- 
arately, and considers how it will move as the wave prop- 
agates. The line traced out by the movement of any one 
of these points is a characteristic. The equation associ- 
ated with eqn. @ that governs the characteristic lines 
of E is 



dE_ 
~dt 



v m {E) 



dE 

dz 



0. 



(3) 



with the velocity v m {E) being that for a point on the 
wave with field strength E, which is given by 



Vm{E) 



c 

no 



l + mx^E^/nl 



-1/2 



(4) 



In this picture, the nonlinearity manifests itself by giving 
different (fixed) velocities to characteristics of different 
E. Thus, as the wave profile E(t) propagates forward in 
space (z), the wave profile becomes distorted by temporal 
compressions or expansions. A shock occurs if a region is 
compressed to the point where two characteristics inter- 
sect. For a nonlinearity, characteristics with higher 
E 2 move more slowly, dragging the peaks toward later 
times; a shock will first occur on the profile where E 2 is 
changing most rapidly in time. Waves approaching the 
point of shocking are shown in fig. [1] 




FIG. 2: Here two points A and B on the field profile (with 
fields Ea, Eb) follow their characteristics as the wave propa- 
gates. Separated initially by a time difference dt, they travel 
at different speeds (v and v — dv), and meet at point C. 

Using eqn. Q along with the construction shown in 
fig. we can derive a simple formula for the distance 
to shocking. The figure shows two characteristics AC 
and BC, originating from points A and B, and converg- 
ing towards a shock at C after a distance of L. We have 
drawn the case where the speed associated with AC (rep- 
resented by its gradient) is lower than that of BC. From 
the geometry of the figure, it is easy to show that 



dv 
~dt 



v 

T 



(•5) 



where t, v, and L = vt are respectively time, velocity, 
and distance. Differentiating the velocity leads to 



dv r , 



d(E m ~ 1 ) 



dt 



2(l + mx {m) E m - l /nl) 
mx {m) v 3 d(E m ~ 1 ) 



3/2 



dt 



2c 2 m dt 
and combining eqns. © and (J7J) yields 



dv m _ 2cn oy /l + mx^E" 1 - 1 /^ 
~dT ~ mx ( - m) (-dE m - 1 /dt) 



(6) 
(7) 

(8) 



Thus L m depends on E(t), and so will vary across the 
pulse profile. As the profile propagates and evolves, so 
will the predictions for L m , nevertheless, the point at 
which shocking occurs relative to the origin remains fixed. 
For a given profile, a shock will occur first at the point 



2 



Kinsler-2007arXiv-mocnl 



MOCNL 



Dr.Paul.Kinsler@physics.org 



where the L m reaches its minimum value. We can there- 
fore define the shocking distance, for any arbitrary wave- 
form E(t), as 



where 



5,, 



2cnQ 



Min 



c m / x (m) 



, (9) 



(-dE m ~ 1 /dt)_ 
\Jl + mx (m) -E m ~V n o- (10) 



We see that the shocking distance is inversely propor- 
tional to the nonlinear strength ^( m ) (as would be ex- 
pected), and that the other important quantity is a 
derivative of powers of the field, i.e. (dE" 1 " 1 /dt). Note 
that for x^" 1 -* > 0, the nonlinear correction factor C m will 
always be greater than 1 for odd-order nonlinearities, but 
may be less than 1 for even-order nonlinearities. In the 
limit of weak nonlinearity, mx' m ' E" 1 ^ 1 jn\ <C 1, 



2cuq 



rn 



Min 



i/x (m) 



(-dE m - l /dt) 



(11) 



III. UNI-DIRECTIONAL MODEL 



the motion of its characteristics. Note also that eqn. 
(H3J) is equivalent to a first-order expansion of the square 
root term in eqn. (j4j. We can again use the method of 
characteristics. Differentiating the velocity u m gives 



du m 
~dt 



no mcx^ m ^ /2n, 



1 



mx 



(m) El 



■V2n§) 



dt 



-(15) 



and combining eqns. 
distance 



Sm+ — 



2cn 



© and (fTS"]) yields the shocking 



Min 



(-dE^/dt) 



(16) 



This is the same as the weak nonlinearity limit of the 
bi-directional prediction in eqn. @, i.e. it lacks the 
correction factor C m defined in eqn. (JTDJ) . 

It is important to note that eqn. (|16p was found us- 
ing an explicitly uni-directional formalism, in which the 
only approximation was of uni-directional propagation 
with no coupling to the backward wave. Thus it allows 
an unambiguous comparison of uni-directional and bi- 
directional propagation. 



Using the directional field approach of Kinsler et al. Q 
we can write a pair of coupled first order wave equations 
for directional fields under the influence of an m- 
th order nonlinearity without dispersion. In the present 
case of a dispersionless medium, these are defined by us- 
ing G ± = JiE x ± y/JIEHy, where e = e (l + x (1) )- This 
combination of scaled fields provides us with a G + field 
whose Poynting vector points forward along the z axis, 
and a G~ field whose Poynting vector points backward. 
To get a uni-directional model, we simply set the back- 
ward propagating G~ field to zero, leaving us with a sin- 
gle first order wave equation for the forward propagating 
G + field. Importantly, we do not require the use of an 
exponential carrier function to impart the directionality 
(see e.g. [13, H3]), as this is achieved by the construction 
of the G ± fields themselves. As in section [Til the calcu- 
lation below can easily be generalised to include a sum 
of nonlinear terms of different order, if so desired. 

After defining E ± = G ± /2- v /e [p, the coupled bi- 
directional wave equations are 

8E ± , n Q dE± = ± x^d(E+ + E_) m 



c dt 



2c 



at 



Setting E- = gives us a single, uncoupled, uni- 
directional wave equation which we can rewrite as 



dE^ 



oz 



0. 



with velocity 



u m (E + ) = — l + mx^E^^nl 



(13) 



(14) 



Comparing eqn. (flU)) to eqn. we see they have 
the same form; so that this wave equation also describes 



IV. DISCUSSION 

We can see by comparing the shocking distances pre- 
dicted by the (exact) bi-directional theory in eqn. @ 
and the (approximate) uni-directional theory in eqn. (|16D 
that both theoretical predictions have the same two dom- 
inant trends: they are inversely dependent on the nonlin- 
ear coefficient x^ m \ an d on the gradient of the (m— l)th 
power of the field. 

The difference between the bi- and uni-directional the- 
ories lies in the nonlinear correction factor C m , which 
only becomes significant for extremely strong nonlinear- 
ities. For small nonlinearities, and correspondingly long 
shocking distances, the poorly phase matched backward 
component does not build up, so the forward field (and 
hence shocking distances) are barely affected. In fused 
silica near the damage threshold, effective nonlinearities 
of order x (3) # 2 ^ 0.06 can be achieved. In the absence 
of dispersion, this amount of nonlinearity would lead to 
shocking distances of less than three wavelengths (e.g. 
S3 ~ 2^m for A = 800nm), a regime in which C3 ~ 1.04. 

For stronger nonlinearities, significant conversion can 
take place before the field has propagated even one wave- 
length, let alone the several needed for the averaging ef- 
fect caused by poor phase matching. The growth of a 
significant backward component then affects the prop- 
agation; indeed it can even be very strongly altered if 
mx^ m ' E" 1 ^ 1 /rip ~ 1. However, such extreme nonlinear- 
ities are not achievable in realistic systems, since they 
require field intensities far beyond the damage thresh- 
olds of standard nonlinear materials. Nevertheless, the 
point of this paper is to examine the limits of the uni- 
directional approximation: adding more realistic mate- 
rial models would test only those models, not the uni- 
directional approximation. 
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FIG. 3: Scaled shocking distances for x' 3 ' as a function of 
nonlinearity: MOC predictions for the exact bi-directional 
case (solid line) uni-directional approximation (dotted line). 
Simulation results are denoted by symbols: with E,H (□), 
G (+), uni-directional G + (*); with distances determined 
using LDD shock detection. 
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FIG. 4: Scaled shocking distances for as a function of 
nonlinearity: MOC predictions for the exact bi-directional 
case (solid line), uni-directional approximation (dotted line). 
Simulation results are denoted by symbols, with E, H (□), 
G (+), uni-directional G + (*); with distances determined 
using LDD shock detection. 



To support the theory, I have also done simulations 
of the bi-directional and uni-directional cases using the 
PSSD technique @. These are (a) straightforward sim- 
ulations of Maxwell's equations, which are naturally bi- 
directional; (b) bi-directional simulations of Maxwell's 
equations, using the G ± directional fields [H; (c) uni- 
directional simulations, using a G + -only forward prop- 
agating model. I treat nonlinear orders m = 2, 3 only 
and propagate two-cycles of a (sinusoidal) CW field sam- 
pled by 512 or 2048 points, giving time resolutions of 
dt = 0.0128fs and 0.0032fs. For the E, H simulations I 
use a Yee-style [22j staggered grid and the PSSD method 
[j| ; for the G ± or G + simulations I use a leapfrog method, 
which, being analogous to the staggered grid, ensures 
that the numerical performance of the simulations is com- 
parable. Numerical "shocks" are detected by using the 
local discontinuity detection (LDD) technique [l7l ]. 

The shocking distances vary slightly from the compa- 
rable theory for two reasons. First, the LDD technique 
is not a perfect predictor of shocks. Second, the stag- 
gered grids used by the simulations make it very hard 
to perfectly match the two initial fields values needed, 
especially in the case of strong nonlinearity. 

Figs. [3] shows the behaviour of the shocking distance 
for the usual case of nonlinearity, and fig. 2] shows it 
for the alternative case of x^ nonlinearity; both show ex- 
cellent agreement between the MOC theory and the simu- 
lations. Note that the bi-directional simulations confirm 
the bi-directional MOC theory, and the uni-directional 
simulations confirm the uni-directional MOC theory. In 
particular, the bi-directional simulations using G ± fields 
agree remarkably well with the Maxwell's equations sim- 
ulations (as they should); this serves to emphasise the dif- 
ference between these and the contrasting uni-directional 
simulations based on (only) the G + field. 

Note that these results are based on analytic results ob- 
tained using the MOC, which cannot account for material 
dispersion. Consequently, while it is still possible to use 
them to validate (or not) the use of a uni-directional ap- 



proximation in a physical model, as I do here, it is hard to 
see how they could be experimentally tested. Any experi- 
ment would instead need to compare its results with sim- 
ulations of the chosen nonlinear material that included 
linear dispersion, an accurate model of the nonlinear re- 
sponse, and perhaps even transverse effects. The test 
would then be a match with a bi-directional simulation, 
but a disagreement with a uni-directional simulation. 

V. THE ROLE OF THE BACKWARD FIELD 

Figs. [3] and |4] exhibit different trends for the scaled 
shocking distances X Ss an d X a difference due to 
the nonlinear correction C m , which contains a field de- 
pendent part F m = mx'" 1 ' , E m ~ 1 / 'rig. In the even-order 
case F m can take either sign; whereas for odd-order it can 
only have the sign of x ■ Thus for sufficiently strong 
nonlinearities, C' m could become zero, implying a zero 
shocking distance. This can happen for odd-order non- 
linearities only if x^ is negative, but will always be pos- 
sible in the even-order case. This exotic "instant shock" 
behaviour is a mathematical prediction rather than a 
physical one, since the effective nonlinearities required 
are far in excess of those attainable in experiment, and 
realistic nonlinearities do not have an instantaneous re- 
sponse. Nevertheless, it is instructive to examine the 
reason for this surprising effect, as it demonstrates how 
backward waves can influence the forward ones. 

We see this instant shock regime in fig. QJ where the 
trend for x &2 is downward, with a sudden dip toward 
zero when the peak field Eq gives F% = — 1, and remain- 
ing at zero for still stronger nonlinearities. This causes 
numerical difficulties, leading to relatively poor agree- 
ment between theory and simulation in this region on 
fig- HI although this improves as the temporal resolution 
is increased. 

On fig. [5] we see the forward and backward field com- 
ponents G ± in the case of a positive x' 3 \ where x^ 3 ^3 
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FIG. 5: Comparison of forward (G + ) and backward (G ) 
field contributions at the LDD shocking point for a x' 3 ' non- 
linearity with C3 ~ 1.32. The G~ field has been scaled up to 
enhance detail. 
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FIG. 6: Comparison of forward (G + ) and backward (G ) 
field contributions at the LDD shocking point for a x' 2 ' non- 
linearity with C2 — 0.71. The G~ field has been scaled up to 
enhance detail. 

increases with nonlinearity. Both are steepest near ~ lfs, 
where G~ has a gradient of the opposite sign to that of 
G + , so the gradient of E ~ G + + G~ is reduced. Con- 
sequently the backward G~ wave has the effect of in- 
creasing the shocking distance; although the effect would 
reverse for a negative x ■ Note that both these compo- 
nents are strongly coupled to each other. 

In contrast, in the x case on fig- G3 the equivalent 
region (near 1.3fs) shows that the gradients of G + and 
G~ can have the same sign, enhancing the gradient of E. 
Closer inspection shows that the gradient of G~ abruptly 
switches sign where G + has a point of inflection, so the 
G _ field mitigates shocking on one side, and enhances it 
on the other. The enhancement then acts like a feedback 
process, where the backward wave steepens the gradi- 
ent, which in turn enhances the backward wave, and so 
on. Thus not only does the backward wave decrease the 



shocking distance, but for sufficiently strong nonlineari- 
ties the feedback can cause a shock in an infinitesmally 
short distance - just as predicted by eqn. (fT6|) . 

Note that since the effect of the x*- 2 -* nonlinearity de- 
pends on the sign of the field, the backward field com- 
ponents will tend to increase the shocking distances in 
some regions of the field profile, but reduce them in oth- 
ers. However, since we always look for the minumum 
shocking distance, this will always be reduced, whatever 
the sign of x^ ■ 



VI. CONCLUSIONS 

I have demonstrated the fundamental limits of the 
widely used uni-directional propagation approximation. 
This was done by comparing analytic results for the 
shocking distance obtained from both an exact bi- 
directional model, and an approximate uni-directional 
model. These theoretical results were for simple nonlin- 
earities of arbitrary order, based on the method of char- 
acteristics, and are supported by numerical simulations of 
both models. The exact bi-directional results were based 
on Maxwell's equations using E & H . The approximate 
uni-directional results relied on the construction of the 
G^ directional fields Q , which enabled the forward-only 
approximation to be made without introducing any ad- 
ditional assumptions. 

The results show that the condition 
Imx^" 1 -* E m ~ x /uq\ <§C 1 must hold for the uni-directional 
approximation to be true; even when no backward field 
is initially present. This condition is usually easily 
satisificd in nonlinear optical materials - even with fields 
strong enough to approach the damage threshold. If 
this condition does not hold, then significant non phase- 
matched forward-backward coupling can occur, affecting 
the propagation accordingly; in extreme cases causing 
behaviour like the "instant shock" discussed in section 
Ivl Such features, demonstrated in a comparison between 
two models which only differ by a uni-directional approx- 
imation made in one, provide an important indication of 
the limitations of a uni-directional description. 
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Appendix: Time propagated MOC 

1. Bi-directional form 

Here we take the expression in eqn. ([3]), and swap the 
roles of t and z. The associated equation governing the 
characteristic lines of E(z) is 



, . dE dE 
w m (E)— + — = 0. 
at oz 

with the (inverse) velocity w m {E) given by 



w m {E) = — 
c 



l + mx ( m )E m - 1 /nl 



1/2 



(17) 



(18) 



Using eqn. (118p along with as time-space swapped ver- 
sion of the construction shown in fig. [21 we can derive a 
simple formula for the time to shocking. From the geom- 
etry, it is easy to show that 



dw 

dz 



w 

z 



w 



(19) 



where z, w, and T = wz are respectively distance, inverse 
velocity, and time. Differentiating the inverse velocity 



leads to 
dw m n mx 



(m) 



dz 2c n 2 



(l+mx {m) E m - l /nt) 



-i/2 dE ri 



dz 



ln§mx (m) ^dE" 1 " 1 



2 c 2 n 2 

mx (m) _i dE" 1 - 1 
-w„ 



dz 



2c 2 ~ m dz 
and combining eqns. (fl~9l) and ([22]) yields 
,dw m 



"(20) 
(21) 
(22) 



W m/- 



dz 

2c 2 / 'dE" 1 - 1 



mx 



(m) 



dz 



(23) 
(24) 



For a given profile, a shock will occur first at the point 
where T reaches its minimum value. We can therefore 
define the shocking time T m , for any arbitrary waveform 

E(z), as 

T m = ^MinSSt, (25) 



where 



mc dE m - 1 /dz 1 
C m = y/l + mx^E" 1 - 1 /^ (26) 
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C m ~ 1 + -mx^E^/nl (27) 

In the spatially propagated form (eqn. ([9])), the correc- 
tion term is C m , not the C m seen here. 

Note that the size of the first order correction term is 
independent of the nonlinear strength, and depends only 
on the properties of the field profile: 



S7„ 



2ng ^ 3mx (m) B m_1 ^B™" 1 ^ 1 



mcx 



(m) 



2rig \ dz 



6ngmx (m) ^ TO_1 /'/£'"' 1 



2n\mc X ^ m ^ \ dz 

m ~ 1 ' 'dE m - ls ~ " 



2. Uni- directional form 



(29) 
(30) 



and combining eqns. ([19]) and (|22|) yields 

, duw+. 



— w m+/~ 



« ! 7, 



2n c /dE r ' 



(37) 
(38) 



For a given profile, a shock will occur first at the point 
where T reaches its minimum value. We can therefore 
define the shocking time T m + , for any arbitrary waveform 
E(z), as 



2 2n c f dE r 



2r? 3 D 2 /v( m ) 



(39) 



where 



D 2 



mr dE^^/dz' ^ 
l + m X (m) E m ~ 1 /2n 2 (41) 
l+m.x (m) £ m "Vng. (42) 



Here we take the expression in eqn. ([3]), and swap the 
roles of t and z. The associated equation governing the 
characteristic lines of E(z) is 



w 



(E) 



dE dE 



= 0. 



dt dz 

with the (inverse) velocity w m +(E) given by 



(E) = — l + m X (m) E m - 1 /2n 2 



(31) 



(32) 



Using eqn. (|18[) along with as time-space swapped ver- 
sion of the construction shown in fig. [21 we can derive a 
simple formula for the time to shocking. From the geom- 
etry, it is easy to show that 



dw 

dz 



w 
z 



w 



(33) 



where z, w, and T = wz are respectively distance, inverse 
velocity, and time. Differentiating the inverse velocity 
leads to 



dw 



np m X 



(m) 



dz 



2n 2 



(l + m X {m) E m - 1 /2n 2 



o dE 



m — l 



dz 



no 



c 2n\ dz 
mx (m) dE m ~ x 



2noc 2 



dz 



-(34) 
(35) 
(36) 



In the spatially propagated form (eqn. ([TCI) ), the D m - 
like term is simply 1. 

Note that the size of the first order correction term is 
independent of the nonlinear strength, and depends only 
on the properties of the field profile: 



57„ 



2nl ^ mx {m) E m-l f dE m- I 



mc X 



(m) 



dz 



2n 1 mx (m) -E" 
n^mc X ^ m ^ 



m—l 



( dE n 



\ dz 



2n a E 



m — l 



f dE 



m — l 



\ dz 



(43) 
(44) 
(45) 



3. Comparison: Bi- vs Uni- directional 

The comparison is not as neat as for the spatially prop- 
agated case, but nevertheless is very similar. For the 
small-nonlincarity limit (i.e. m X ( m ) E" 1 ' 1 /uq <C 1), the 
difference is just the factor of 1 + m X ^ m ^> E" 1 ^ 1 /2n 2 that 
appears for the spatially propagated results in the same 
limit. 

Note that the shocking time T is not just related to 
the shocking distance L by a simple factor of c, since in 
the spatially propagated picure, the peak-dragging effect 
corresponds to a time offset which needs to be added to 
the propagation time L/c. 
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